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Abstract 

A classical optics waveguide structure is proposed to simulate resonances of short range 
one-dimensional potentials in quantum mechanics. The analogy is based on the well known 
resemblance between the guided and radiation modes of a waveguide with the bound and 
scattering states of a quantum well. As resonances are scattering states that spend some 
time in the zone of influence of the scatterer, we associate them with the leaky modes of a 
waveguide, the latter characterized by suffering attenuation in the direction of propagation 
but increasing exponentially in the transverse directions. The resemblance is complete 
since resonances (leaky modes) can be interpreted as bound states (guided modes) with 
definite lifetime (longitudinal shift). As an immediate application we calculate the leaky 
modes (resonances) associated with a dielectric homogeneous slab (square well potential) 
and show that these modes are attenuated as they propagate. 


1 Introduction 

In general, a resonance may be defined as the excitation of a system that results when one 
of its characteristic frequencies matches a particular value that is defined by either the 
boundary conditions, the external forces, or any other interaction or constraint applied to 
the system. The concept arises from the study of oscillating systems in classical mechanics 
and extends its applications to physical theories like electromagnetism, optics, acoustics, 
and quantum mechanics, among others p[J. In quantum and nuclear physics the resonances 
are associated with metastable states of a system which has sufficient energy to break up 
into two or more subsystems [2j. The whole range of scattering experiments provides 
a big number of examples in this matter m- This last is because scattering includes 
the situation in which the impinging particle is ‘captured’ for a while in the zone of 
influence of the scatterer, so that the projectile and the scatterer form a new unstable 
system (resonance state). The capturing occurs for specific energies of the projectile 
that are defined by the general properties of the scatterer. Eventually, the projectile is 
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released and escapes from the interaction zone (decay process). In this model the ‘time of 
capture’ corresponds to the lifetime of the decaying system that is formed of the scatterer 
and the projectile under the resonance condition. Collisions, in addition, are modeled 
as interactions localized in time and space. This implies that the involved potential 
vanishes rapidly enough in space, so that incoming and outgoing asymptotic states can be 
represented by wave packets in free motion. The problem is usually reduced to the analysis 
of one-dimensional effective potentials such that the main information is obtained from 
the transmission and reflection amplitudes of the scattering states, see e.g. |5j. In the case 
of one-dimensional square potential wells and barriers the resonances appear for specific 
values of the parameters that define the interaction. Indeed, for such parameters the 
resonances can be associated with the bell-shaped peaks of the corresponding transmission 
coefficient W- The presence of resonances can be also stimulated either by shallowing a 
square well [8] or by adding to it a static electric held [9]. As resonances are special cases of 
scattering states, they are represented by irregular (not normalizable) vectors [TOUTT] that 
satisfy the purely outgoing condition [T2] (see also m- This last property of resonances 
is useful in constructing complex potentials that operate as optical devices with position- 
dependent complex refractive index [6] [T3] ITT]. 

Remarkably, quantum wells and barriers can be made with high precision and com¬ 
positional control by modern epitaxial growth techniques [15] . This has opened a wide 
window of experimental possibilities to explore various quantum phenomena in real de¬ 
vices and actual laboratory settings. For example, it has been shown that the trans¬ 
mission time of a wave packet that is tunnelling through a one-dimensional barrier does 
not depend on the barrier thickness pa]; a similar property is found in one-dimensional 
potential wells [13 dg. This condition implies anomalously large group velocities of the 
wave packets that impinge the barrier (well) since the velocity of the packet’s peak must 
increase with the width across the interaction zone. Such a superluminal phenomenon has 
been observed using electromagnetic analogues for evanescent modes [19] and microwave 
pulses [20], at the same time this result has stimulated the designing of high-speed de¬ 
vices based on the tunneling properties of semiconductors (see, e.g., Chs. 11 and 12 of 
Ref. [21]). On the other hand, as the scattering by quantum wells attenuates the outgoing 
wave packets only because of the multiple reflections at the well boundaries, it is expected 
that the resonance condition implies a maximum in the time spend by the projectile in 
traversing the zone of influence of the scatterer [6]. Recent theoretical research includes 
the case of a rectangular well that is embedded in an environment formed of a zero po¬ 
tential energy (flat potential) at the right and a parabolic potential at the left of the 
well [221 - 124] . Thus, the connection between time delay and ‘time of capture’ associated 
to the scattering processes would be verified in the laboratory by using semiconductor 
materials. 

Another accessible laboratory tool to study quantum phenomena is provided by waveg¬ 
uide structures. The connection is based on the analogies between light waves and 
Schrodinger wave-state functions that have been studied since the early years of quantum 
physics [25] • I 11 this approach the amplitude of an electromagnetic held is associated to the 
probability amplitude of a quantum state with two degrees of freedom while the propaga- 
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tion of such field through the waveguide corresponds to the time-evolution of the quantum 
state. Such optical-quantum analogy is representative of a wider scheme that includes 
diverse connections between classical and quantum systems [26], and has been successfully 
applied in studying Anderson localization m, coherent enhancement [2E| and destruction 
of tunneling [2D], coherent population transfer [2D], and decay of metastable states cou¬ 
pled to semi-infinite, tight-binding lattices [213, among others. Of special interest for the 
purposes of this paper, nanoseconds delay-time measurements have been performed into 
the electromagnetic framework by using microwave setups of narrowed waveguides [32], 
The relevance of this last result obeys the fact that the waveguide simulates a quantum 
potential barrier so that the measured electromagnetic delay-time is in correspondence 
with the quantum tunneling time. Thus, the waveguide structures represent an emergent 
and versatile arena to perform realistic experiments addressed to get deeper insights in 
the understanding of resonances in quantum physics. 

The aim of this paper is to get a classical optics analogy of the resonances associ¬ 
ated to one-dimensional short range potentials in quantum physics. The centerpiece of 
our approach is the connection between the mathematical structure of the propagation 
of electromagnetic signals in waveguides and the dynamics of a time-dependent quantum 
state, lu this form our interest is twofold: we want to connect the theoretical calculations 
of resonances with mensurable properties of an electromagnetic wave that propagates in 
a waveguide, and we would like to extend the application of the resonance mathematical 
techniques to the solving of problems associated to electromagnetic waveguides. A col¬ 
lateral result will be the connection between the resonance properties of electrons that 
are injected into heterostructures and the properties of the electromagnetic signals that 
propagate in optical waveguides. The paper is organized as follows. In Section [2] we 
revisit the main properties of the resonance states associated to one-dimensional short 
range potentials. In particular, it is shown the relevance of the Fock-Breit-Wigner distri¬ 
bution in a scattering process. When the transition amplitude from an initial scattering 
state to its time-evolved version is ruled by such a distribution one immediately identifies 
energy resonances between the projectile and the scatterer. It is also shown that analyti¬ 
cal continuation of the eigenvalue equation to a complex plane of energies is required to 
get an appropriate mathematical description of resonance states. In Section [3] we include 
some of the well known generalities of the mode fields in waveguide structures and their 
connection with the Schrodinger equation of a system of two degrees of freedom. The case 
of electromagnetic propagation through a finite, dielectric homogeneous material and its 
association with a one-dimensional square well potential is considered in Section 13.11 as 
immediate application. Section [4] is devoted to the connection between resonance states 
and a singular form of modal field in waveguides, the latter is known as leaky mode and 
is characterized by suffering attenuation in the direction of propagation but increasing 
exponentially in the transverse directions. We apply the resonance approach of Ref. [6] in 
the calculation of the leaky modes associated to the dielectric homogeneous slab described 
in Section 13.11 As an example we provide the eigenvalues (resonances) corresponding to 
the leaky modes supported by a slab of definite width and refractive index, and show 
that these modes are attenuated as they propagate. Section [5] includes the calculation 
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of the longitudinal shifts for the leaky modes of Section [4j we show that the maxima of 
these shifts are in correspondence with the maxima of the transmission coefficient, a re¬ 
sult that verifies the connection between times of capture and resonances in the scattering 
processes. Finally, Section |H] includes some of the main conclusions of the present work. 


2 One-dimensional scattering states and resonances 

Let H be the Hamiltonian defined by a short range one-dimensional potential U(x). 
Assuming that the spectrum of H is formed of discrete and continuous eigenvalues one 
can consider the following sets of state-vectors 


( \4>k), k £ z D c z + 

discrete eigenvalues E k < I D = Yj^x-d |0fc)(0fc| (!) 

( H D = Span{|0 fe )} fcG x c 

( | 4>e), Eei E cm 

continuous eigenvalues A •<->•< I E = Iegx e ^E\ ( I > e){4 > e\ (2) 

l He = Span{|^ E )} E6 x B 

where Ip and I E are respectively the identity operator in Ho and He- Consistently, the 
spectrum and Hilbert space of H are cr(H) = {E k } ke z D U {E} e &i e and H = Ho ®H E 
respectively. As usual, the discrete and continuous eigenvectors ([I])-© are orthonormal 
in their respective subspace according to the rules 

= 8jk, (MM) = HE - E'). (3) 

Hereafter we pay detailed attention to the scattering states (i.e., vectors in He)- The 
description of bound states (i.e., vectors in Hd) is straightforward by using the well 
known relationship f dE -H- Yk- The reason is that, as we have said in the introduction, 
resonance states arise quite naturally in the scattering processes. 

Any scattering state | ijj) G He can be expressed in the basis of continuous eigenvectors 


= 1* 


dEC(E)\M, C(E) = (M1>), 


' EgI e 


and is such that 

\m 2 = m)= [ dE\c( e )\ 2 . 

J E£X e 

In ^-representation the above expressions read as 


'ip(x) = / dEC(E)(j) E (x), C(E) = / dxcj) E (x)ip(x), 

EdzI-E J M. 


(4) 

(5) 

( 6 ) 
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and 

M* = [ dx\j>(x)\ 2 = [ dE\C(E)\ 2 , (7) 

J R J E(zXe 

with ~z standing for the complex conjugate of z G C. The above expressions allow to write 
the energy distribution c o(E) as follows 

w(£) = = |C(£)| 2 . ( 8 ) 

Let (0]) 7 equivalently (ED, be the initial state of a stationary system. The transition 
amplitude T(t > 0) from | ijj(t = 0)) = |Vh) to = | ip t ) is given by the inner product 


T{t > 0) = = / dxij) o(x)^t(x) = / dEu(E)e lEt/h . 


(9) 


' EgXe 


It is clear that T can be investigated in terms of spatial coordinates x or as a function 
of the energy distribution (jHJ). In this work we shall assume that c o(E) is given by the 
following expression 


u(E) 


_(W_ 

(E - E„ + iT/ 2 ) (E-E„- ir/2) 


(r/2 ) 2 

(E - E 0 f + (r/2) 2 ’ 


( 10 ) 


where E 0 G Ie is the energy of the initial state j^o) an d T > 0 is a parameter having 
units of energy. Note that the bell-shaped function (fTOlh it is centered at E = E 0 , for 
which oj(E 0 ) = 1, and has a half-width at half-maximum equal to T/2. In addition, oj(E) 
is a meromorphic function that has two poles E± = EqEiT /2. This function is known as 
either Cauchy (mathematics), Lorentz (statistical physics) or Fock-Breit-Wigner -FBW- 
(nuclear and particle physics) distribution. Using (JS]) we can write the Fourier coefhcient 
C(E) in the expansion (|4]) as follows 


C(E) 


r/2 

E — Eq T iT j 2 


( 11 ) 


After the analytic continuation of Z# C M to the complex U-plane one can introduce (TTUT) 
into © to arrive, up to a global constant factor, at the expression 


T(t > 0) 


-ie 0 t/h 

2 


^ e -iE 0 t/h^j e ~rt/(2h)^ 



( 12 ) 


This last result has the form of a transient oscillation [T] and means that the transition 
from |-0 O ) t° |V’t) is an exponential decreasing function of time whenever C(E) is given by 
Eq. (1IB . That is, the initial state j^o) decays according to the time-dependent probability 

IUf>0)| 2 = (Tj e- r >/«, (13) 

As we can see, the width F of the distribution (fTOl) represents an indirect measure of 
the lifetime of |^o) because at time r := k/Y the probability that this state has not yet 
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decayed is reduced approximately in 36%, i.e., Y 2 /4 —$■ T 2 /(4e). It is clear that Y and r 
are correlated, the smaller the value of Y the larger the lifetime r. In general, one says 
that there is a resonance Eq of width Y when either the energy distribution of the system 
at time t > 0 is given by the FBW function (HUD or the Fourier coefficients C(E) of the 
expansion (JUJ) are given by (HID . 

We would like to emphasize that the dependence of oj, C and T on the complex 
number e 0 = E 0 — iY/2, as this has been expressed in Eqs. (HUD - (HUD . is not merely 
aesthetic; this corresponds to the formulation of Gamow, Fock, Breit and Wigner in the 
late 1920s and along the 1930 decade (see details in HD- It is due to Siegert jT2j that 
a resonance is introduced as a solution u £ of the Schrodinger equation associated to the 
complex eigenvalue e = E — iY j 2 and satisfying the purely outgoing condition [6j: 

lim /3 = =|=i/c, (14) 

#—>•±00 

where k 2 = 2 me/h 2 G C, and 

/3:=--y-ln u £ . (15) 

ax 

It can be shown that the above conditions lead automatically to the FBW distribution 
m with k a number in the fourth quadrant of the complex plane C, see [Tj. In turn, 
the distribution of the complex points e 1 ' 2 = hk/y/2m characterizes a Riemman surface 
of two sheets, each one cutting along the positive real axis. These complex numbers 
define directly the behavior of the propagator (fTUjh For instance, consider a very narrow 
distribution ui(E), that is Y « 1 (the point £ G C is close to the positive real axis). As 
T defines the vicinity E G (E 0 — Y/2,E 0 + Y/2), we see that ‘large lifetimes’ r = h/Y 
means narrow widths T, so that the energy E of the system at time t > 0 is very close 
to the initial energy E 0 . In the extremal case Y —y 0, the distribution c o(E) is as narrow 
and singular as S(E — E 0 ). That is, the related lifetime r is large enough (r —> +oo) to 
consider the initial vector |^ 0 ) a non-decaying (stable) state. In this context, bound states 
are considered as decaying states with infinitely large lifetimes and vice versa, unstable 
(decaying) states correspond to bound states with a definite lifetime. 

To get more insights about the complex eigenvalue £ 0 let us consider the Schrodinger 
equation 

H*l>(x;t) — ih—ij}{x\t). (16) 

Assuming separation of variables ip{x\t) = <p(x)Q(t ), with Q(t) = e - lEt l n an d E G Xe the 
separation constant, the probability density does not depend on time |^(x;t)| 2 = |<^(a;)| 2 . 
Now, let us make the analytic continuation of the eigenvalue equation H<p(x) = E<p(x) in 
order to include complex eigenvalues £ = E — iY /2, that is Hu e = eu £ . In this case we have 
t) = u £ (x)e~ l£t / h : and the probability density is an exponential decreasing function of 
time t)| 2 = \u £ (x)\ 2 e~ rt ^ h . Clearly, the spatial part u e (x) of the wave-function ^(x; t ) 
is no longer a square-integrable function because this belongs to the complex eigenvalue 
£ and the Hamiltonian H is still Hermitian. Indeed, it can be proved that u e (x) diverges 
exponentially as |x| —>■ +oo. Thus, the density |'0(x;t)| 2 increases exponentially for either 
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large |x| or large negative values of t. The usual form to avoid some of the complications 
connected with the limit t —» — oo is to consider the long lifetime limit T —» 0, see e.g. [33]. 
Therefore, we consider the condition 


T/2 

AE 


« 1 . 


(17) 


This last means that the level width T must be much smaller than the level spacing AE 
in such a way that closer resonances imply narrower widths (longer lifetimes). In general, 
the main difficulty is precisely to find the adequate E and T. However, it has been 
shown [6] that for one-dimensional stationary short range potentials the superposition 
of a denumerable set of FBW distributions (each one centered at each resonance E n , 
n — 1,2...) entails an approximation of the transmission coefficient T such that the larger 
the number N of resonances involved, the higher the precision of the approximation: 


N 

T ph u j n (E) = ^u;(E,E n ); 

71 = 1 


u(E,E n ) 


(IW2) 2 

(E - E n ) 2 + (T n /2) 2 ' 


(18) 


3 Mode fields in waveguide structures 


Consider the Helmholtz equation 

^2 

—E(x, z) + -rp^E(xi z ) + kln 2 (x)E(x, z ) = 0, (19) 

with k 0 = w/c the wave number in vacuum. Equation (TT9l) can be achieved from the 
Maxwell equations for transversal electric fields E in presence of inhomogeneous dielec¬ 
tric materials, neglecting the magnetic polarizability and the dispersive properties of the 
medium. Indeed, (1IU1) corresponds to an electric wave propagating along the positive 
^-direction and polarized in the ^-direction, assuming that the refractive index n depends 
only on the x-coordinate. In the paraxial regime we use the ansatz 

E(x, z) = ip(x)e~ iko ^- no)z , (20) 


to transform (jT9|) into the linear second order differential equation 


1 d 2 
2 A)q77. 0 dx 2 



<p(x) = £<P (x), 


£ = 1 ~n 0 , 


( 21 ) 


known as the paraxial Helmholtz equation. In (120]) - ([2T]) the number hq is a reference 
refractive index, that may be taken as the maximum of n(x), and we have assumed that 
the weakly guiding condition (An = |n — no| << 1) holds. In this approach the parameter 
e, known as the propagation constant, is the value of the ^-component of the generalized 
linear momentum [34] and defines the slope of the light beam as follows 


e = —p z = —n{x) cos 6{x). 


( 22 ) 
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In general, modal fields associated to open waveguides are classified into guided 
(trapped) modes and radiation modes, the latter propagating outside the core [351438] . 
The guided modes are bounded fields that correspond to real values of the propagation 
constant fulfilling — n 0 < e < —1. They are perfectly guided through the waveguide as 
they suffer no attenuation along the optical axis and decay exponentially in the transverse 
direction. Thus, as guided modes are trapped into the core of the waveguide they add 
no energy flux into the infinite surrounding clad. In turn, propagating radiation modes 
correspond to real propagation constants fulfilling —1 < £ < 0 , so that they are bounded 
oscillatory fields in the transverse as well as in the longitudinal directions (35]. They 
appear as a consequence of the infinite extension of the clad cross-section. Guided and 
radiation modes constitute a complete set of modal fields in the sense that any physical 
solution to the paraxial Helmholtz equation can be expressed as a linear superposition of 
a finite number of guided modes and a continuum of radiation modes. 

At this stage it is important to stress the resemblance of the electric held expression 
(1201) with the solutions t ) = ip(x)e~ lEt i h of the Schrodinger equation (fT 6 |) . This allows 
the identification 

k^ 1 oh, zPl, e = 7 — no O E. (23) 

In adittion, Eq. (j2H would correspond to the eigenvalue equation Hip = Eip, provided 
the following additional identification holds 

fk~ 2 \ d 2 

n 0 o m, —nix) Hl/(i), H — — f -2— ) —— — nix). (24) 

\2tt-o ) dx 2 

The identification (l23li - (j24li is well known in the literature, see e.g. [34] • This allows to 
construct the classical optics analogy of diverse quantum phenomena and vice versa, the 
quantum analogy of classical optical phenomena [261134] . The analogy includes the sets 
of vectors introduced in (JT|) and ([ 2 ]) because bound states correspond to ‘trapped’ waves 
that have no chance of escaping from the zone of influence of the scatterer (as we have 
seen in the previous section, these states can be interpreted as decaying states with in¬ 
finitely large lifetime). Thus, quantum bound states correspond to guided electromagnetic 
modes. Particles in scattering states, on the other hand, can be temporally trapped by the 
scatterer in such a form that their presence in the zone of interaction decreases in time. 
In this form, scattering states correspond to radiation modes. The orthonormalization 
properties of guided and radiation modal fields can be stated in a form similar to that 
of the quantum mechanical bound and scattering states (see, e.g. |35]|36]). Hence, the 
total radiated power can be expressed as the sum of two terms: a first one associated to 
the power that propagates inside the core, given by a hnite sum containing the contribu¬ 
tion of guided modes, and a second term corresponding to the power radiated from the 
core, expressed as an integral over all radiation modes. In this form, denoting by A *. and 
A(e) the guided and the radiation mode amplitudes, respectively, the total power of the 
radiated held in the waveguide axis direction is given by |37) : 



On the other hand, as we have seen, resonances are scattering states of complex eigen- 
energy that spend definite intervals of time in the interaction zone. So resonances are not 
connected to radiation modes because these last do not belong to complex eigenvalues. 
Next we shall take full advantage of the analogy indicated above to get optical models for 
the resonances discussed in Section [21 


3.1 Wave propagation through a dielectric homogeneous slab 


As an immediate application let us consider the simplest case of electromagnetic propa¬ 
gation through a finite, dielectric homogeneous material. It is convenient to rewrite the 
eigenvalue equation (T2Th as follows 


d 2 


k ~ 2 


2 noJ dx 2 


+ 1 — n(x ) 


ip= (e + l)ip- 


( 26 ) 


Considering a combination of media distributed along the x-axis such that there is vacuum 
in (—oo, — a) U (a,+oo), and a medium of constant refractive index n — Uq > 1 in 
x G [—a, a], the equation (126]) is in correspondence with the eigenvalue equation defined 
by a linear square well of depth Vo = |1 — Uq\ and width 2a. The solutions to this equation 
are well known and can be found elsewhere, thought the analysis of resonances and the 
resolution of the transmission coefficient in terms of the FBW distribution, see Eq. 
can be consulted in Ref. (6j. Hereafter we shall label the regions of the media (vacuum, 
medium, vacuum) from left to right as /, II and III. 


4 Leaky modes and resonances 

Besides the guided and radiation modes discussed in the previous section, it is well known 
the existence of solutions to the Helmholtz equation with complex propagation constants 
that have an “unusual behavior”: they increase exponentially in the transverse direction 
and suffer attenuation in the longitudinal one [35]. These solutions are called leaky modes 
and have attracted much attention over the time. The leaky modes do not belong to 
the set of guided and radiation modes of infinite transverse waveguides, though they 
constitute a discrete set of modes that can be used to approximate the continuous mode 
contribution to the radiation held, in regions close to the core and far from the source, 
in order to avoid integration over all radiation modes [3&H40] . Yet, leaky modes have 
properties that are very similar to the ones of guided modes, but the former modes loss 
energy as they propagate along the optical axis [36]. In slab waveguides this leakage is 
due to refraction at the waveguide core-clad boundaries (refracting leaky modes [37]). 
In waveguides with curved cross-section the leakage can be also produced by the partial 
reflection at the boundaries due only to the curvature of the cross-section (tunneling 
leaky modes) ]35jl37]. It can be shown that the leakage is slower in the tunneling modes 
compared to the refractive modes [35]. In any case, these modes behave much like guided 
modes since they are able to propagate long distances inside the guide before their power 
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loss started to be significant. Different methods have been applied to determine the 
leaky modes of planar [36][4T| as well as cylindric [35t[39j waveguides involving isotropic 
and anisotropic dielectrics [42]. These techniques include numerical methods as mode¬ 
matching [43], finite differences and finite element [3511421144] . approximate methods as 
WKB [45], transfer matrix [4T] and differential transfer matrix methods [401146] . There is a 
conceptual discussion concerning leaky modes, the exponential growing in the transversal 
direction of these nearly guided modes make them appear as nonphysical solutions to the 
Helmholtz equation. Yet, this behavior is required if the energy flux must be preserved: 
As the flux is damped in the longitudinal direction, there must be an increasing flux in 
the transverse direction [36]. 


The problem of finding leaky modes is mathematically equivalent to determine the 
resonant states in quantum mechanics for a particle in a short range interaction potential. 
This is convenient because we may apply different methods that have been addressed to 
obtain the resonance energies and the corresponding wave functions in order to determine 
the complex propagation constants and the modal fields. To this end, let us regard 
the square well potential analogy indicated in Section 13.11 with a complex propagation 
constant, e = £r — i The corresponding complex wave number k = kn + ikj is related 
to £ through 

k'f{ ~ k'j T k n kj 

£R ~ 2 kl ’ 2 ~~ kl 

The energy density flux in the transverse direction is not null because k r ^ 0 (see [ 6 ]). A 
wave coming from region III enters into resonance with the optical medium and is thus 
trapped inside the core (region //), then it propagates and travels a finite distance along 
the optical axis. After a while, the wave is finally emitted into the clad. The result is 
a wave propagating to the left in region /, and another one propagating to the right in 
region III. The correct sign for kji is thus positive and the whole process obey the pure 
outgoing condition. The expression for the electric energy density fulfills 



lim p e oc e k ° rz \ip(x)\ 2 , (28) 

x —>-±oo 


meaning that the power decreases exponentially as 2 : grows whenever T > 0 and ki < 0. 
The optical-quantum analogy requires also a purely outgoing wave <p(x), see Section [2] so 
the amplitudes of the external electric fields increase exponentially as x —> ± 00 . Leaky 
modes then correspond to wave numbers k = kjj + ikj in the fourth quadrant of the 
complex plane. The factor k^ 1 /T can be interpreted as the mean distance traveled by the 
beam, along the optical axis, before its intensity decays by 1/e of its intensity at a focus. 
In this context T represents a measure of the transient confinement of the field into the 
core. In the limit T —)■ 0 we will have a quasi-stationary mode propagating inside the 
slab in a very similar way as a guided mode. In this limit, it is possible to establish an 
analytic approach to calculate the propagation constant. Bearing this in mind, following 
CO}. we may write 


r/2 

A£_r 


<c 1, 


(29) 
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2(fco) Z(k 0 ) Z(ko) 



Figure 1: Propagation of resonances inside a slab waveguide of width k^a = 30 and constant refractive 
index n = Uq = 1.5. The slab supports 17 leaky modes (24 < m < 40) defined by Equations (15(1 and 
m- The figure shows the propagation in the axial direction (from the bottom to the top in each slab) of 
the real part of the leaky modes field amplitudes for m = 24 (left), to = 32 (center) and to = 40 (right). 
Note that these are attenuated as the mode propagates. The leakage is faster for larger values of to. 


where is the separation between two consecutive resonances. To get the position £r 
and the width T of the resonances we follow the approximations reported in [E]. After 
some calculations one arrives at the expressions 


£r 


1 (™Y~u 0 


and 


2/cq(7o V 2 a J 
r i 


2 k 0 aU 0 

with m an integer number fulfilling 


\/2 (£r + 1 ), 


\/2U 0 (U 0 - 1) < 




, < J2U2. 

2k 0 a v 0 


(30) 

(31) 


(32) 


Thus, for each combination of the parameters (a,Uo,ko), there is a finite number of 
resonances given by (130T) and (I3TT) . Table [0 shows the eigenvalues e corresponding to the 
17 leaky modes supported by a slab waveguide of refractive index Uq = 1.5 and width 
k 0 a = 30. The behavior of the real part of the field amplitude of these modes, as they 
propagate in the longitudinal direction through the waveguide, is depicted in Figure Q] for 
m = 24 (left), m = 32 (center) and m = 40 (right). Observe that the modes with smaller 
values of m attenuate slower compared to those with larger values of m. 


5 Longitudinal shifts 

Let us determine the longitudinal shift of a wave packet under scattering by an optical 
medium. This will provide us some insight of the distance traveled by the leaky mode in 
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k 0 a = 30, 

Uo 

= 1.5 


m 

£ = £r 

,k 

1 2 

m 

£ = £r 

L 2 

24 

-0.973621 - 

zO.0051042 

33 

-0.504815 - 

zO.0221150 

25 

-0.928842 - 

zO.0083833 

34 

-0.443587 - 

zO.0234424 

26 

-0.882236 - 

zO.0107847 

35 

-0.380531 - 

zO.0247350 

27 

-0.833802 - 

zO.0128120 

36 

-0.315647- 

zO.0259981 

28 

-0.783540 - 

zO.0146215 

37 

-0.248936 - 

zO.0272358 

29 

-0.731450 - 

zO.0162860 

38 

-0.180397- 

zO.0284514 

30 

-0.677533 - 

zO.0178462 

39 

-0.110031 - 

zO.0296476 

31 

-0.621788 - 

zO.0193273 

40 

-0.037836 - 

zO.0308267 

32 

-0.564215 - 

zO.0207462 





Table 1: Eigenvalues e corresponding to the 17 leaky modes supported by a slab of width k^a = 30 and 
constant refractive index Uo = 1.5. 


the waveguide before the longitudinal power is completely lost. This problem is equivalent 
to evaluate the reflection and transmission phase times for a particle in a rectangular well 
potential pH [T7J HI]. Consider a wave-packet scattering solution for the external electric 
field 

E I (x,z)= I dkf(k)e i{ - kx ~ ko£{ - k)z) + 

Jo 

E ni (x,z)= [ dfc/(fc)f(fc)e^ fea: - 6 W feoz - 2fea+ ^^ +7r/2 ), (34) 

Jo 

where f{k) is the Fourier coefficient, k is defined by e(k) = 0, r(k) and t[k) are the 
refiexion and transmission amplitudes respectively and 

<j){k) = - arctan ^ = U ° ^ + 2k ° (U ° ~ X )) ' ( 35 ) 


dkf(k)r( /t) e*( _fc3;_£ ( fc ) fc o 2 - 2 fca+ A fc )) 


(33) 


Applying the stationary phase condition [ISl[23|[24] to the incident and transmitted wave 
packets one has 


d_ 

dk 


k 2 

kx - —z + koz 

2 k 0 


0, 


d_ 

dk 


k 2 Tj- 

kx — 2 ka - —z + k 0 z + </>(/c) H— 

2k 0 2 


0. (36) 


So that the equations of motion for the centers of the corresponding Fourier components 
are given by 


x = 


k 

kt 


-z, x = 2 a + 


k dcp 

ko dk 


(37) 


The point Zi n at which an incident plane wave of wave number k would reach the left wall 


of the optical medium is then Zi n = — ^-a. In the same way the point Zt at which the 
reflected or transmitted wave leaves the optical medium is 




( 38 ) 
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Sz/a 5z/a 




(a) (b) 

Figure 2: Longitudinal shift as a function of the eigenvalue e = —p z for a dielectric slab of width koa 
equal to (a) 30 and (b) 50,000. In both cases Uo = 1.5. 


Then the longitudinal shift of the reflected or transmitted wave is readily written 


Sz = zt — Zi, 


k 0 dp 
k dk 


(39) 


Sz/a T 



Figure 3: Longitudinal shift (continuous line) and transmission coefficient (dotted line) of an electric 
wave propagating in a dielectric slab of width koa = 30, and constant refractive index Uo = 1.5, both as 
a function of the eigenvalue e. Observe that the maxima of both curves coincide. 

Figure [2] shows the longitudinal shift of a Fourier component under reflection or transmis¬ 
sion as a function of £ for (a) k 0 a = 30 and (b) k^a = 50, 000. It is worthwhile to point 
out that the peaks of this shift coincide with the peaks of the corresponding transmission 
coefficient as shown in Figure [3] This implies that the electromagnetic modes having max¬ 
imum transmission probabilities, and thus corresponding to the resonant (leaky) modes, 
travel a maximum longitudinal distance in the core before they are finally transmitted to 
the infinite clad. Figure 0] shows the longitudinal shift as a function of the width k 0 a for 
the eigenvalue £ = —0.995 and Uq = 1.5. We can appreciate that the transmission shift 
Sz takes negative values, these have been already predicted for both rectangular barriers 


13 



































5 10 15 

kftCL 


Figure 4: (a) Longitudinal shift as a function of the width k^a of a dielectric slab with constant refractive 
index Uq = 1.5. Here e = —0.995. 


and wells in the quantum approach mum. The analogy between the Schrodinger and the 
paraxial Helmholtz equations suggest that this phenomenon may be used to measure the 
negative phase time in the scattering of particles by a potential well (compare with [T8]). 


6 Conclusions 

The classical optics analogy of quantum systems is based on the relationship between 
the Helmholtz equation and the two degrees of freedom Schrodinger equation. Although 
such analogy is mathematical, the optical solutions are useful in simulating the quantum 
ones for diverse cases. For instance, it is well known that the guided and radiation modes 
of a waveguide structure can be put in correspondence with the bound and scattering 
states of a one-dimensional short range potential in quantum mechanics. However, the 
association of leaky modes with their quantum counterpart, as well as the identification 
of the optical analogue of a quantum resonance, is rarely discussed in the literature on the 
matter. One of the reasons would be that in both cases the calculation of the eigenvalues 
and wave-solutions is not as immediate as in the scattering (radiation) or bound (guided) 
states (modes). In this work we have presented a form to connect these two kinds of 
wave solutions: a leaky mode is in correspondence with a resonance state. The former 
is as close to a guided mode as a resonance is close to a bound state. That is, a leaky 
mode (resonance state) is a guided mode (bound state) for which the longitudinal shift 
(lifetime) is definite. 
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